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Abstract 



We study the evolution of electromagnetic fields in curved spacetimes and calculate the 
exact wave equations of the associated electric and magnetic components. Our analysis ap- 
plies to a general spacetime, is fully covariant and isolates all the sources that affect the 
CO ■ propagation of these waves. Among others, we explicitly show how the different parts of 

the gravitational field act as driving sources of electromagnetic disturbances. When applied 
to perturbed FRW cosmologies, our results argue for a superadiabatic-type amplification of 
large-scale cosmological magnetic fields in Friedmann models with open spatial curvature. 

O PACS numbers: 04.20.-q, 98.80.-k, 41.20.Jb 

o 

"o ! 1 Introduction 

$h . Electromagnetic studies in curved spaces have long established the direct coupling between the 

Maxwell and the Einstein fields. The interaction emerges from the vector nature of the electro- 
magnetic field and from the geometrical approach to gravity introduced by general relativity and 
it is interpreted as a sort of scattering of the electromagnetic waves by the spacetime curvature. 

In the present article we study electromagnetic fields in general curved spacetimes by using 
the covariant approach to general relativity. Our analysis is non-perturbative, in the sense that 
it does not perturb away from a given metric but provides the full nonlinear equations before 
linearising them about a chosen background. In addition, we study the physically measurable 
electric and magnetic components of the Maxwell field, rather than the Faraday tensor or the 
electromagnetic 4-potential. This on one hand complements earlier work on the subject, while 
on the other it allows for a more compact mathematical presentation and for a more transparent 
physical interpretation of the results. The evolution of the Maxwell field is studied in a general 
spacetime without imposing any a priori symmetries on the latter. The only restriction is on 
the matter component which is of the perfect fluid form. We derive, from first principle, the 
electric and magnetic wave equations and identify all the kinematical, dynamical and geometrical 
sources that drive the propagation of these waves. We demonstrate the effects of the observers' 
motion and show how the different parts of the gravitational field, namely the Ricci and Weyl 
fields, affect propagating electromagnetic disturbances. Moreover, the non-perturbative nature 
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of our formalism means that the nonlinear equations apply to variety of situations of either 
astrophysical or cosmological interest. 

With the full equations at hand, we proceed to consider the evolution of electromagnetic fields 
in spatially curved Friedmann-Robertson- Walker (FRW) models. Noting that the symmetries of 
the FRW spacetime are generally incompatible with the generic anisotropy of the electromagnetic 
field, we consider the evolution of the latter in perturbed Friedmann universes. In particular, 
we look at the spacetime curvature effects on the linear evolution of the magnetic component 
of the Maxwell field. Our results show that, when the model is spatially closed, the magnetic 
field has an oscillatory behaviour with a decreasing amplitude according to the familiar aT 2 - 
law (where a is the FRW scale-factor). The same is also true for spatially open models with 
the exception of large-scale magnetic fields. There, we find that the field decays as a -1 and 
therefore that magnetic flux conservation no longer holds at long wavelengths. This result 
can be seen as an effective superadiabatic-type amplification of large-scale magnetic fields in 
spatially open FRW universes due to curvature effects alone. Crucially, the amplification is 
achieved without introducing any new physics and without breaking away from the standard 
properties of Maxwell's theory. 

We start with an outline of the covariant formalism in section 2 and provide a covariant 
treatment of the electromagnetic and gravitational fields in sections 3 and 4 respectively. The 
nonlinear electromagnetic wave equations are derived in section 5 and in section 6 they are 
linearised and solved around curved FRW models. We discuss our results in section 7. 



2 The covariant description 

The covariant approach to general relativity uses the kinematic quantities of the fluid, its energy 
density and pressure and the gravito-electromagnetic tensors instead of the metric which in itself 
does not provide a covariant description. The key equations are the Ricci and Bianchi identities, 
applied to the fluid 4-velocity vector, while Einstein's equations are incorporated via algebraic 
relations between the Ricci and the energy-momentum tensors. Here, we will only give a brief 
description of the approach and direct the reader to a number of review articles for further 
details and references [l]-[4]. 



2.1 The 1+3 spacetime splitting 

Consider a general spacetime with a Lorentzian metric g a b of signature (— , +, +, +). Then, 
allow for a family of fundamental observers living along a timelike congruence of worldlines 
tangent to the 4-velocity vector 

dx a 

Ua = "dT ' ^ ' 

where r is the associated proper time and u a u a = — 1 [4]. This fundamental velocity field 
introduces an local, 1+3 "threading" of the spacetime into time and space. The vector u a 
determines the time direction and the tensor h a b = g a b + u a Ub projects orthogonal to u a into 
what is known as the observers' instantaneous rest space. Note that, in the absence of rotation, 
h a b also acts as the metric of the spatial sections. 
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Employing u a and h ab one defines the covariant time derivative and the orthogonally pro- 
jected gradient of any given tensor field T a b... cd "' according to 

•/;,/, .."'- = u e V e T ab .. cd - and B e T ab .. cd - = h e s hj h^h/ ■ ■ ■ V s T fp .J r - , (2) 

respectively. The former indicates differentiation along the timelike direction and the latter 
operates on the observers' rest space. 

2.2 The matter field 

Relative to the aforementioned fundamental observers, the energy-momentum tensor of a general 
imperfect fluid decomposes into its irreducible parts as [4] 1 

Tab = l^u a u b + ph ab + 2q^ a u b) + n ab . (3) 

Here, fi = T ab u a u b and p = T ab h ab /3 are respectively the energy density and the isotropic 
pressure of the medium, q a = —h a b T bc u c is the energy-flux vector relative to u a and ir ab = 
h( a c hb) d T C( i is the symmetric and trace- free tensor that describes the anisotropic pressure of 
the fluid. 2 It follows that q a u a = = ir ab u a . When the fluid is perfect both q a and 7r ab are 
identically zero and the remaining degrees of freedom are determined by the equation of state. 
For a barotropic medium the latter reduces to p = p(fi), with c 2 = dp/dp representing the 
associated adiabatic sound speed. 

When dealing with a multi-component medium, one needs to account for the velocity "tilt" 
between the various matter components and the fundamental observers [5]. Here, however, we 
will consider a single-component fluid and we will assume that the fundamental observers are 
moving along with it. 

2.3 The covariant kinematics 

The observers' motion is characterised by the irreducible kinematical quantities of the u a - 
congruence, which emerge from the following covariant decomposition of the 4-velocity gradient 

VbU a = <Jab+Uab + \®h ab -U a Ub, (4) 

where a a b = D^ a )> ^ab = ^[b u a]^ © = V a ii a = D a u a and ii a = u b V b u a are respectively the shear 
and the vorticity tensors, the expansion (or contraction) scalar and the 4-acceleration vector [4]. 
Then, a ab u a = = u) a bU a = u a u a by definition. Also, on using the orthogonally projected 
alternating tensor e a b c (with e a b c = 3u^ a e bc ^ d u d ) , one defines the vorticity vector to a = e abc uj bc /2. 

The nonlinear covariant kinematics are determined by a set of three propagation equations 
complemented by an equal number of constraints [4]. The former contains Raychaudhuri's 
formula 

6 = -|6 2 - §(p + 3p) - 2{a 2 - uj 2 ) + B a u a + u a u a , (5) 

1 Throughout this article we use geometrised units with c = 1 = 8irG. Consequently, all geometrical variables 
have physical dimensions that are integer powers of length. 

2 Angled brackets denote the symmetric and trace-free part of projected second-rank tensors and the orthogo- 
nally projected component of vectors. 
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for the time evolution of 0, the shear propagation equation 

<j{ab) = -§@0"a6 ~ Vc{a° C b) ~ ^(a^b) + D( a U b ) + U( a U b ) ~ E ab + \lT ab , (6) 

which describes kinematical anisotropics, and the evolution equation of the vorticity 

Co( a ) = -§6w a - §curh2 + a ab uj b . (7) 

Note that a 2 = a ab a ab /2 and oj 2 = uj ab u> ab /2 = oj a uj a are respectively the magnitudes of the 
shear and the vorticity, while E ab is the electric component of the Weyl tensor (see section 3.2 
below). Also, cxa\v a = e abc D b v c for any orthogonally projected vector v a by definition. 

Equations (5), (6) and (7) are complemented by a set of three nonlinear constraints. These 
are the shear 

DVafc = §D a 6 + curlu; a + 2e abc u b uj c - q a , (8) 

the vorticity 

B a uj a = ii a uj a , (9) 

and the magnetic Weyl constraint 

H ab = cwc\<T ab + T) {a uj b) + 2ii {a uj b) , (10) 
where curlT a 5 = e cd ^ a D c T b } d for any orthogonally projected tensor T ab . 

3 The electromagnetic field 

Covariant studies of electromagnetic fields date back to the work of Ehlers [1] and Ellis [3]. In 
addition to its inherent mathematical compactness and clarity, the covariant formalism facilitates 
a physically intuitive fluid description of the Maxwell field. In particular, the latter is represented 
as an imperfect fluid with properties specified by its electric and magnetic components. 

3.1 The electric and magnetic components 

The Maxwell field is covariantly characterised by the antisymmetric electromagnetic (Faraday) 
tensor F ab , which relative to a fundamental observer decomposes into an electric and a magnetic 
component as [3, 6] 

F ab = 2u [a E b] + e abc H c . (11) 

In the above E a = F ab u b and H a = e abc F bc /2 are respectively the electric and magnetic fields 
experienced by the observer. Note that E a u a = = H a u a , ensuring that both E a and H a 
are spacelike vectors living in the observer's 3-dimensional rest-space. Also, expression H a = 
e abc F bc /2 guarantees that H a is the dual of the antisymmetric (pseudo) tensor F ab . 

The Faraday tensor also determines the energy-momentum tensor of the Maxwell field. In 
particular we have 

T^ b m) = -F ac F c b -\F cd F cd g ab , (12) 
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which, on using (11), provides an irreducible decomposition for . More precisely, relative 
to a fundamental observer, the latter splits into [3, 6] 

T ir } = U E2 + H 2 )u a u b + \{E 2 + H 2 )h ab + 2Q (a u b) + V ab . (13) 

Here E 2 = E a E a and H 2 = H a H a are the magnitudes of the two fields, Q a = e abc E b H c is the 
electromagnetic Poynting vector and V a b is a symmetric, trace-free tensor given by 

V ab = V(ab) = \(E 2 + H 2 )h ab - E a E b - H a H b . (14) 

Expression (13) provides a fluid description of the Maxwell field and manifests its generically 
anisotropic nature. In particular, the electromagnetic field corresponds to an imperfect fluid 
with energy density [E 2 + H 2 )/2, isotropic pressure (E 2 + H 2 )/6, anisotropic stresses given by 
V a b and an energy-flux vector represented by Q a . Equation (13) also ensures that ri em) a = 0, 111 
agreement with the trace- free nature of the radiation stress-energy tensor. Finally, we note that 
by putting the isotropic and anisotropic pressure together one arrives at the familiar Maxwell 
tensor, which assumes the covariant form 

M ab = \{E 2 + H 2 )h ab - E a E b - H a H b . (15) 
3.2 Maxwell's equations 

We follow the evolution of the electromagnetic field by means of Maxwell's equations. In their 
standard tensor form the latter read 

V [c F ab] =0, and V b F ab = J a , (16) 

where (16a) manifests the existence of a 4-potential and J a is the 4-current that sources the 
electromagnetic field, With respect to the u a -congruence, the 4-current splits into its irreducible 
parts according to 

Ja = PeU a + J a , (17) 

with p e = —J a u a representing the charge density and J a = h a b J b the orthogonally projected 
current (i.e. J a u a = 0). 

Relative to a fundamental observer, each one of Maxwell's equations decomposes into a 
timelike and a spacelike component. Thus, by projecting (16a) and (16b) along and orthogonal 
to the 4- velocity vector u a , we obtain a set of two propagation equations [3, 6] 

E {a ) = {°ab + e abc io c - l9h ab ) E b + e abc ii b H c + curlF a - J a , (18) 

H {a ) = ((Tab + e abc uj c - l@h ab ) H b - e abc u b E c - curLE a , (19) 
and the following pair of constraints 

B a E a + 2uj a H a = Pe , (20) 

B a H a - 2co a E a = . (21) 

Note that in addition to the usual "curl" and "divergence" terms, there are terms due to the 
observer's motion. According to Eq. (20), in the absence of an electric field the observed charge 
density is p c = 2uj a H a - This means nonzero charge density unless uj a H a = (see [7] for a 
discussion on the charge asymmetry of the universe). Also, following (21), the magnetic vector 
is not a solenoidal unless to a E a = 0. 
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3.3 The conservation laws 

The antisymmetry of the Faraday tensor (see Eq. (11)) and the second of Maxwell's formulae 
(see Eq. (16b)) imply the conservation law 

V a J a = 0, (22) 

for the 4-current density. Then, on using decomposition (17), expression (22) provides the 
covariant form of the charge density conservation law [3, 8] 

p c = -6pe - D a J a - u a J a . (23) 

Thus, in the absence of spatial currents, the charge density evolution depends entirely on the 
average volume expansion (or contraction) of the fluid element. 

3.4 Ohm's law 

The electrical conductivity of the medium determines the relation between the 4-current and 
the associated electric field via Ohm's law. In covariant form the latter reads 

Ja ~ PeU a = oE a , (24) 

where a is the scalar conductivity of the medium [9]. Projecting the above into the observer's 
rest space one arrives at 

J = oE a . (25) 

Thus, non-zero spatial currents are compatible with a vanishing electric field as long as the 
conductivity of the medium is infinite (i.e. for a — > oo). Alternatively, one can say that at the 
infinite conductivity limit, which defines the well known MHD approximation, the electric field 
vanishes in the frame of the fluid. On the other hand, zero electrical conductivity implies that 
the spatial currents vanish even when the electric field is non-zero. 

4 The gravitational field 

Covariantly, the local gravitational field is monitored by a set of algebraic relations between the 
Ricci curvature tensor and the energy- momentum tensor of the matter. The free gravitational 
field, on the other hand, is described by the electric and magnetic components of the conformal 
curvature (Weyl) tensor. 

4.1 The local Ricci curvature 

In the general relativistic geometrical interpretation of gravity, matter determines the spacetime 
curvature which in turn dictates the motion of the matter. This interaction is manifested in the 
Einstein field equations, which in the absence of a cosmological constant take the form 

Rab = T ab - \Tg ab , (26) 
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where R ab = R ac b c is the spacetime Ricci tensor, T ab is the energy- momentum tensor of the 
matter fields, with T = T a a being the trace. For our purposes the total energy-momentum 
tensor has the form T ab = + T^ m \ where is given by Eq. (3) and T^ rn) by Eq. (13). 
Thus, 

T ab =ln+±(H 2 + E 2 )] u a u b +[p + \(H 2 + E 2 )} h ab + 2{q (a + Q {a )u h) + ir ab + V ab , (27) 

ensuring that \i + (H 2 + E 2 )/2 is the total energy density of the system, p + (H 2 + E 2 )/6 is 
the total isotropic pressure, q a + Q a is the total heat flux vector and ir ab + V a b is the total 
anisotropic pressure. The inclusion of electromagnetic terms in the energy-momentum tensor of 
the matter guarantees that the contribution of the Maxwell field on the spacetime geometry is 
fully accounted for. 

Starting from the Einstein field equations and assuming that T ab is given by Eq. (27), we 
arrive at the following algebraic relations [8] 

R ab u a u b = ±(p + 3p + E 2 + H 2 ), (28) 
h a b R bc u c = -(Qa + Qa), (29) 
h a c h d R cd = [\(n-p+\E 2 + \H 2 )]h ab + TT ab + V ab . (30) 

In addition, the trace of (26) gives R = — T, with R = R a a and T = T a a = 3p—fi, where the latter 
result is guaranteed by the trace- free nature of T^™ 1 ^ . Note that the above expressions are valid 
irrespective of the strength of the electromagnetic components. When the Maxwell field is weak 
relative to the matter, namely for E 2 , H 2 <C /x, one might treat the electromagnetic contribution 
to the spacetime curvature as a first order perturbation. Finally, recall that q a = = 7r ab when 
dealing with a perfect fluid. 



4.2 The long-range Weyl curvature 

The Ricci tensor describes the local gravitational field of the nearby matter. The long-range 
gravitational field, namely gravitational waves and tidal forces, propagates through the Weyl 
conformal curvature tensor. The splitting of the gravitational field into its local and non-local 
components is demonstrated in the following decomposition of the Riemann tensor 

Rabcd = Cubed + \ (gacRbd + QbdRac ~ QbcRad ~ 9adRbc) ~ \R {.9ac9bd ~ 9ad9bc) , (31) 

where C abc d is the Weyl tensor. The latter shares all the symmetries of the Riemann tensor 
and is also trace-free (i.e. C c acb = 0). Relative to the fundamental observers, the Weyl tensor 
decomposes further into its irreducible parts according to 

C a bcd = {9abqp9cdsr ~ VabqpVcdsr) U^U E P — {j] a bqp9cdsr + 9abqpVcdsr) U^U H P , (32) 

where g abcd = g a c9bd ~ 9ad9bc (e.g. see [10, 11]). The symmetric and trace-free tensors E ab and 
H ab are known as the electric and magnetic Weyl components and they are given by 

E a b = C ac bdU c u d , and H ab = \e a cd C c d be u e , (33) 
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with E a bU = = H a bU . Given that E a b has a Newtonian counterpart, the electric part of the 
Weyl tensor is associated with the tidal field. The magnetic component, on the other hand, has 
no Newtonian analogue and therefore is primarily associated with gravitational waves [2]. Of 
course, both tensors are required if gravitational waves are to exist. For a comparison with the 
non-covariant metric based treatments of gravitational waves we note that in perturbed FRW 
models the harmonically decomposed, pure-tensor metric perturbation is Ht = 2E + cr'/n [16]. 
Here, E and a represent the harmonic parts of the transverse traceless electric Weyl and shear 
tensors respectively. Also, n is the associated wavenumber and a prime denotes derivatives with 
respect to conformal time. 

The Weyl tensor represents the part of the curvature that is not determined locally by matter. 
However, the dynamics of the Weyl field is not entirely arbitrary because the Riemann tensor 
satisfies the Bianchi identities. When contracted the latter take the form [10] 

V'Cabcd = V[ fe i? a ] c + i g c [bV a ]R , (34) 

by means of decomposition (31). In a sense the contracted Bianchi identities act as the field 
equations for the Weyl tensor, determining the part of the spacetime curvature that depends on 
the matter distribution at other points [10]. The form of the contracted Bianchi identities guar- 
antees that once the electromagnetic contribution to the Ricci curvature has been incorporated, 
through the Einstein field equations, the effect of the Maxwell field on the Weyl curvature has 
also been fully accounted for. 

Expression (34) splits into a set of two propagation and two constraint equations, which 
monitor the evolution of the electric and magnetic Weyl components [2]- [4]. These formulae 
are not used to derive the electromagnetic wave equations of section 5.3 and are therefore not 
essential for our purposes. Here we simply note that the aforementioned set of equations is 
remarkably similar to Maxwell's formulae, which in turn explains the names of E a b and H a b. 
This Maxwell-like form of the free gravitational field underlines the rich correspondence between 
electromagnetism and general relativity, which has been the subject of theoretical debate for 
many decades (see [12]-[15] for a representative list). 

5 The electromagnetic wave equations 

Studies of electromagnetic waves in curved spacetimes have long established that, while propa- 
gating similar to any other travelling wave, electromagnetic disturbances also interact with the 
spacetime curvature. As a result, electromagnetic signals propagate inside as well as on the fu- 
ture light cone of an event, indicating the failure of Huygens' principle in curved spaces [17]-[19]. 

5.1 The wave equation for the electromagnetic field tensor 

Maxwell's equations immediately provide a wave equation for the electromagnetic field tensor. 
In particular, taking the covariant derivative of (16a) and using (16b) we arrive at 

V 2 F ab = -2R acbd F cd + R a c F cb + F a c R ch + V b J a - V a J b , (35) 

where V 2 = V a V a is the generalised covariant Laplacian operator (e.g. see [20, 21]). The above, 
which holds in a general spacetime, reveals the role of the curvature as a driving source of 
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electromagnetic disturbances. Note that the Riemann and Ricci curvature terms in the left- 
hand side of Eq. (35) emerge after using the Ricci identity 

2V[ V b ]-Pcd = R a bceF e d + RabdeFc , (36) 

which here monitors the commutation between the covariant derivatives of F a b. Expression (35) 
can also provide the individual wave equations for the electric and magnetic components of F^. 
For example, contracting Eq. (35) along u a eventually leads to the wave equation of E a , while its 
dual provides the magnetic wave equation. Here, we will follow an alternative route and obtain 
these expressions directly from the decomposed Maxwell formulae (18) and (19). 

5.2 The electro/magneto-curvature coupling 

In addition to the Einstein field equations, vector sources, like the electromagnetic field, obey 
an extra set of equations, known as the Ricci identities, which manifest the direct interaction 
between electromagnetism and spacetime geometry. This coupling emerges naturally from the 
vector nature of the Maxwell field and from the geometrical approach to gravity of general 
relativity. When applied to the magnetic field vector the Ricci identity reads 

2V [c V b] H a = R dabc H d ; (37) 

with an exactly analogous expression for the electric component. Clearly, on using decomposition 
(31), the Ricci identity couples the electromagnetic field explicitly with both the local and the 
long-range gravitational field. Also, by projecting the above into the observer's rest space one 
arrives at what is known as the 3-Ricci identity 

2D[ c Db]-ff~a = -2e c bdU d H( a ) + TZdabcH d , (38) 

describing the interaction between the magnetic field and the local spatial geometry [22, 23]. 
Clearly an exactly analogous relation holds for E a as well. Note that 1Z a b c d is the orthogonally 
projected part of R a bcd-> namely the Riemann tensor of the observer's local 3-space. Note that 
the validity of both (37) and (38) extends to any arbitrary spacetime (e.g. see [2, 10]). 

5.3 The wave equations for the electric and magnetic fields 

Equations (18)- (19) monitor the propagation of electromagnetic fields in a general spacetime 
either in vacuum (i.e. for source-free fields with p c = = J a ) or in the presence of matter. 
Starting form these formulae one can work out the wave equations for propagating electromag- 
netic radiation in a general spacetime. In particular, taking the time derivative of Eq. (18) one 
obtains the wave-like evolution equation of the electric field. Similarly, the time derivative of 
Eq. (19) leads to the corresponding wave equation of the magnetic field. In the Minkowski space 
of special relativity these calculations are relatively straightforward since the geometry of the 
space is trivial. In the context of general relativity, however, this is no longer true and one has 
to account for the coupling between the electro-magnetic fields and the spacetime geometry dis- 
cussed earlier. Technically speaking, this requires using the Ricci identities and leads to spatial 
curvature terms every time the projected derivatives of E a or H a commute. In addition, the 
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Ricci identities guarantee a Weyl field contribution whenever a time derivative and a projected 
gradient of either the electric or the magnetic field commute. 

Assuming that the matter component has a perfect fluid form with a barotropic equation 
of state, we take the time derivative of Eq. (18) and project it orthogonal to u a . Then, using 
the kinematical propagation and constraint equations of section 2.3, expression (14), relations 
(19)-(32) and the commutation laws (37), (38) we arrive at the following wave equation for the 
electric field vector 

E {a) - D 2 E a = i/x(l + 3w)E a + (a ab - e abc uo c - §0/t a6 ) E b + ±6 (a ab + e abc u; c - f 6h a6 ) E b 
-a {a c a b)c E b + e abc E b a cd co d + § (a 2 - |w 2 ) E a + \u {a ui b) E b + ii b u b E a 
-le abc u b cnv\E c + D {a E b) u b + |e abc F b D c G + e abc H d D b a cd + D {a uj b) H b 
+%e abc H b cm\u; c + 2D {a H b) J - 2e abc a b d D {c H d ^ + e abc u b H c + lu b u b H a 
+lH b uj b u a - 3u b H b oo a + 3e abc u b a cd H d + \ Pc ii a - D a/ 9 C - &J a - j a 
-K ab E b - E ab E b + H ab H b . (39) 

Similarly, one may start from Eq. (19) and proceed in an analogous way to obtain the wave 
equation of the magnetic field vector 

H {a) - B 2 H a = |/x(l + 3w)H a + (a ab - e abc co c - §6/i ab ) H b + \Q (a ab + e abc co c - ±®h ab ) H b 
-a {a c a b)c H b + e abc H b a cd Lo d + § (a 2 - |w 2 ) H a + \u {a u b) H b + ii b u b H a 
-le abc u b CMi\H c + B {a H b) u b - le abc E b D c Q - e abc E d V> b a cd - B {a u b) E b 
-|e a6cJ B 6 curlu; c - 2V {a E b) J + 2e abc a b d T>^E d ^ - e abc u b E c - \ii b u b E a 
-\E b u b ii a + 3ii b E b u a - 3e abc ii b a cd E d - |p e w a + 2e abc ii b J c + curlj" a 
-K ab H b - E ab H b - H ab E b . (40) 

As expected, when there are no charges and currents, one recovers Eq. (40) from (39) by simply 
replacing E a with H a and H a with — E a . Similarly, we obtain (39) from (40) after substituting 
H a with E a and E a with — H a . In the presence of charges and currents, however, this symmetry 
no longer holds and the apparent breakdown reflects the absence of magnetic monopoles. 

The above expressions provide a covariant description of propagating electromagnetic waves 
in a general spacetime and incorporate the electromagnetic input to the curvature of the latter. 3 . 
So far the only restrictions are those imposed on the fluid, which has a barotropic equation of 
state. That aside, Eqs. (39) and (40) are fully nonlinear in perturbative terms. Once the back- 
ground is specified, these equations can describe the evolution of the electromagnetic field at any 
perturbative level. In general, of course, one needs to couple these formulae with the appropri- 
ate propagation equations of the various kinematical, dynamical and geometrical variables that 

By including the Maxwell field in the Einstein field equations (see Eqs. (26)-(30)) the electromagnetic contri- 
bution to the spacetime geometry has been fully accounted for. In practise this means ensuring that /i has been 
replaced with n + (E 2 + H 2 )/2, p with p + (E 2 + H 2 )/6, q a with Q a and -K a b with V a b in every formula used to 
derive Eqs. (39) and (40). For example, by implementing the aforementioned substitution into the kinematical 
expressions of section 2.3, we incorporate fully the electromagnetic impact on the model's kinematics. 



10 



appear in the right-hand side of (39) and (40). Clearly, the more complicated the background 
the more equations are necessary for the system to close. 

Among others, the above given wave equations show how the kinematical quantities, namely 
the expansion, the shear, the vorticity and the acceleration drive the propagation of electromag- 
netic waves. Here, the barotropic nature of the matter component means that the 4-acceleration 
takes the form 

H(l + w)u a = -B aP + p c E a + e abc J b H c , (41) 

with contributions from gradients in the fluid pressure and from the electromagnetic Lorentz 
force only. The input from the spacetime geometry to Eqs. (39) and (40) is through the spa- 
tial and the Weyl curvature components. The former is represented by lZ ab , the orthogonally 
projected 3-Ricci tensor, defined by 

Tlab = K c acb = h a c h b d R cd + R achd u c u d + v ac v c b - ev ab , (42) 

where v ab = T) b u a is the second fundamental form describing the extrinsic curvature of the space 
(e.g. see [10, 23]). Note that the tidal part of the Weyl field contributes to the evolution of either 
E a or H a via its direct coupling with the aforementioned fields. The effect of the magnetic Weyl 
tensor, on the other hand, is indirect and requires the presence of both the electromagnetic field 
components. 

The non-perturbative nature of our analysis, namely the fact that we have not yet specified 
our background spacetime, means that Eqs. (39) and (40) apply to a range of physical situa- 
tions (e.g. see [24]- [26]). For example, in the absence of matter sources one can always set the 
observer's acceleration to zero (see Eq. (41)). If, in addition, the spacetime is stationary and 
non-rotating (i.e. set = = u a ), expression (39) reduces to 

E a ~ V 2 E a = a ab E b - a {a c a b)c E b + \a 2 E a + e abc H d V b a cd - 2e abc a b d B^ c H d > 

-K ab E b - E ab E b + H ab H b , (43) 

with an exactly analogous wave equation for H a . When the shear and the Weyl components 
are divergence- free (i.e. for T) b a ab = = D b E ab = T) b H ab ), the above describes the propagation 
of electromagnetic radiation in the presence of gravitational waves alone. Thus, using Eqs. (43) 
one can look form a different perspective into the age old problem of the interaction between 
electromagnetic and gravitational waves in isolated astrophysical environments, away from the 
gravitational field of massive compact stars (e.g. see [27]- [30] and references therein). In what 
follows, however, we will consider a cosmological application of (39) and (40). 

6 Electromagnetic fields in curved FRW models 

The generic anisotropy of the electromagnetic energy-momentum tensor makes the Maxwell 
field incompatible with the high symmetry of the FRW spacetime. The implication is that 
the simplest models where one can study cosmological electromagnetic fields are the perturbed 
Friedmann universes. 
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6.1 The linear wave equations 

Consider a FRW background cosmology with curved spatial sections. In covariant terms, the 
isotropy of the FRW model translates into uj a = = a ab = u a and E ab = = H ab , while 
their spatial homogeneity ensures that all orthogonally projected gradients vanish identically 
(i.e. D a /i = = D a p = D a 0). This means that fj,, p, 6, lZ ab = lZh ab /3 and their time derivatives 
are the only non vanishing background quantities. 

When studying cosmological electromagnetic fields there is a widespread perception that, 
given the conformal invariance of the Maxwell field and the conformal flatness of the FRW 
spacetimes, flat spaces provide an adequate background (e.g. see [31, 32]). This is only approx- 
imately true however, since the FRW symmetries are generally incompatible with the presence 
of electric or magnetic fields. As it is clearly stated in [33], adopting the conformal triviality 
of Maxwell's equations on FRW backgrounds means ignoring the electromagnetic impact on 
the FRW symmetries. This is a good approximation when dealing with weak electromagnetic 
fields, but only on small scales in models with nontrivial spatial geometry. In the latter case, 
the approximation becomes progressively less accurate as one moves on to larger scales and the 
3-curvature effects start kicking in. Putting it in another way, with the exception of fully in- 
coherent radiation, one must study cosmological electromagnetic fields in perturbed Friedmann 
universes. The latter, however, are no longer conformally flat. 

On these grounds, we consider a perturbed Friedmann universe with non-Euclidean spatial 
sections and allow for the presence a weak electromagnetic field. The latter vanishes in the 
background, thus guaranteeing that both the electric and the magnetic field vectors are first- 
order, gauge-invariant perturbations [34]. Then, the source-free components of the nonlinear 
wave equations (39) and (40) linearise to 

E a - B 2 E a = -l@E a - lG 2 E a + i/i(l + 3w)E a - K ab E b , (44) 

and 

H a - B 2 H a = -\OH a - l@ 2 H a + 1^(1 + 3w)H a - K ab H b , (45) 

respectively. During linearisation quantities with nonzero background value have zero pertur- 
bative order, while those that vanish in the background are first order perturbations and higher 
order terms are neglected. For example, the Weyl-free nature of the FRW metric guarantees that 
the Weyl effects are nonlinear. The 3-Ricci curvature, on the other hand, contributes to both 
(44) and (45). Recall that TZ ab = (2k/a 2 )h ab to zero order, where k = 0, ±1 is the curvature 
index and a represents the scale-factor of the unperturbed model. In other words, the symme- 
tries of the FRW metrics ensure that, to linear order, the electromagnetic field interacts only 
with the 3-Ricci part of the spacetime curvature. The curvature terms in (44) and (45) reflect 
the earlier mentioned coupling between electromagnetism and spacetime geometry. Unless the 
background model is spatially flat, these are clearly first-order perturbative terms and should 
be taken into account in any complete linear study of cosmological electromagnetic fields. These 
linear curvature terms clearly show why large-scale electromagnetic fields are not adequately 
treated on flat FRW backgrounds. 

Given that the source-free E a and H a fields satisfy identical linear wave equations, we will 
only consider the magnetic component and proceed by introducing the following harmonic de- 
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composition for H a 

H a = Y,H (n) Q^\ (46) 

n 

where n is the comoving eigenvalue of the n-th harmonic component and are the asso- 
ciated vector harmonics. As usual D a H^ = = Qa^ and are eigenfunctions of the 
Laplace-Beltrami operator so that D 2 Q^ = — (n 2 / a^Q^. Employing decomposition (46) and 
introducing the conformal time variable rj (with fj = 1/a) we recast Eq. (45) as 

H U + ^(n) = "4 (£) H{ n) - 2 (£) 2 H {n) - 2 (^j H {11) - 2kH (n) , (47) 

where a prime indicates differentiation with respect to rj. Then, on introducing the "magnetic 
flux" variable H( n ) = a 2 H^, the above reduces to 

W(n) + n2 K(n) = -2feW (n) . (48) 

This is a wave equation for H( n ) with a driving term on the right-hand side which depends on 
the background spatial curvature and vanishes only when the background is spatially flat. Note 
that in a model with closed spatial sections the Laplacian eigenvalue is given by n 2 = v(y + 1), 
where v takes the discrete values v = 1, 2, 3, . . . . Alternatively, n 2 = v 2 + 1 when k = — 1 and 
n 2 = v 2 for k = (with z/ 2 > in both cases). 



6.2 The linear solutions 

The driving term in the right-hand side of Eq. (48) is clearly sensitive to the sign of the back- 
ground spatial curvature. Let us consider first a FRW model with closed spatial sections. When 
k = +1, Eq. (48) takes the form 

H'( y) + [2 + v(v+l)]H iv) = 0, (49) 

with v = 1,2,3, ... . The above leads to the following oscillatory solution for the v-th magnetic 
mode 

H {v) = {Ci cos v / 2~T^(VTl>7 + C 2 sin ^2 + v(v + 1)77] } , (50) 

where Ci and C2 are constants. In other words, for k = +1, the magneto-curvature term in the 
right-hand side of (48) does not have any significant effect on the evolution of the field, which 
oscillates in time with an amplitude that decays according to the a _2 -law. The only difference 
relative to the k = case is a change in the oscillation frequency near the long wavelength limit. 
Note that the oscillatory behaviour of the field is ensured on all scales by the compactness of 
the closed space. 

When dealing with the hyperbolic geometry of the spatially open FRW model, however, the 
oscillatory behaviour of H( n ) is not always guaranteed. Indeed, for k = — 1 Eq. (48) takes the 
form 

n'( u) + (u 2 -i)n (v) = o, (5i) 
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with v 2 > 0. Clearly, when v 2 > 1 the harmonic mode T~L( V ) oscillates just like in a perturbed 
closed FRW model. On these scales the background geometry makes no real difference in the 
evolution of the field. This agrees with our perception that curvature effects become progressively 
less important as we move towards smaller scales. On sufficiently long wavelengths (i.e. for 
v 2 < 1), the geometrical effects take over and Eq. (51) no longer accepts an oscillatory solution. 
In particular, as v 2 — > we have 

n {u} = d cosh r] + C 2 sinh 77 = C 3 e v + C 4 e~ v , (52) 

where C\ and C2 are constants and £3,4 = (C± ±£2) /2. Note that, since n 2 = v 2 + 1 > 1 always, 
these long wavelength solutions still correspond to subcurvature modes [35]. To have a closer 
look at the effect of geometry on the linear evolution of the field, we note that the evolution of 
a spatially open FRW model is monitored by 

a© = 3coth(/3r ? ), (53) 

with (3 = (1 + 3ttj)/2 by definition and (3rj > (e.g. see [36]). The above holds throughout 
the various periods in the lifetime of an open FRW universe, provided the barotropic index 
w remains constant during each epoch. Then, the relation between the scale factor and the 
conformal time variable is 

a = a »{\—^) e ™> (54) 

where 770, ao depend on the initial conditions. Throughout the dust era w = and (3 = 1/2, 
while w = 1/3 and (3 = 1 when radiation dominates. Finally, during a period of inflationary 
expansion with p = —p we have (3 = — 1. Note that in the latter case the conformal time variable 
takes negative values. According to expression (54), there are extensive periods in the lifetime 
of the universe (i.e. as long as 77 or 77 » 0) when the relation between the cosmological scale 
factor and the conformal time variable is (see also [37]) 

a oc e v . (55) 

Substituting this result into the right-hand side of Eq. (52), and taking into account that T~t(v) = 
a 2 H^ by definition, we arrive at 

H {v) = da' 1 + C 4 a- 3 . (56) 

Therefore, large-scale magnetic fields in perturbed spatially open FRW models decay 
rate considerably slower than the standard "adiabatic" a _2 -law. The immediate consequence is 
that, at the long wavelength limit, the cosmological magnetic flux is no longer conserved. Instead, 
the product a 2 H^ increases with time. This opens the possibility of an effective superadiabatic 
amplification of the field on large scales similar to that found in [38] . Even if the universe is only 
marginally open today, this effect could have important implications for the present strength 
of primordial large-scale magnetic fields. Particularly for those fields that survived an epoch 
of inflation, since they would be much stronger than previously anticipated. Note that during 
inflation the conductivity of the cosmic medium is effectively zero, which in turn ensures the 
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absence of spatial currents (see Section 3.4). In this article, we have focused primarily on the 
mathematics of the magneto-geometrical interaction and provided a qualitative measure of its 
implications for large-scale magnetic fields. A discussion of the physics, together with a detailed 
quantitative study of the amplification effect, will be given in a subsequent article. 

So far, similar modifications in the evolution of cosmological magnetic fields have been ob- 
tained at the expense of standard electromagnetic properties, and in particular of the conformal 
invariance of the Maxwell's equations (e.g. see [38]- [43] for a representative, though incomplete, 
list). Moreover, in some cases this effect is achieved by introducing new physics. Our analysis 
shows that one can still arrive at the same result by taking into account the natural, general 
relativistic coupling between the electromagnetic field and the spacetime curvature. In other 
words, contrary to the widespread perception, superadiabatic magnetic amplification is possible 
within conventional electromagnetic theory. Here, this has been done through the field's cou- 
pling to the intrinsic curvature of spatially open FRW models. Interestingly, however, analogous 
effects can also occur in perturbed flat FRW cosmologies by coupling the magnetic field to the 
Weyl curvature of the model, namely to the gravitational waves [44] . All these cast new light 
on the role and the potential implications of spacetime geometry for the evolution of large-scale 
cosmic magnetic fields. 

7 Discussion 

The general relativistic coupling between the electromagnetic and the gravitational fields has 
long been known in the literature. So far, this interaction has been primarily studied in terms 
of the Faraday tensor and of the electromagnetic 4-potential [17]-[21]. Here, we have taken an 
alternative approach by looking at the evolution of the individual electromagnetic field com- 
ponents in a general curved spacetime. Assuming that the matter field is of the perfect fluid 
form, we have derived from first principle the nonlinear wave equations of the electric and the 
magnetic parts of the Maxwell field. This complements earlier studies which have provided a 
differential/integral formulation of Maxwell's formulae in terms of the physically measurable 
components of the electromagnetic field (e.g. see [45]- [47]). Our approach identifies and iso- 
lates all the sources that drive the propagation of electromagnetic fields by keeping the separate 
aspects of the problem quite distinct. Also, by being manifestly covariant at every step, our 
calculation avoids undue complexity without introducing any specific coordinate frame. We 
show explicitly how the electric and magnetic fields are affected by the various kinematical and 
dynamical quantities and particularly by the different parts of the gravitational field. 

Given that large-scale electromagnetic fields are generally incompatible with the FRW sym- 
metries, we consider perturbed models and concentrate on the evolution of large-scale magnetic 
fields. In particular, we linearise our equations about spatially curved FRW spacetimes and 
investigate the implications of the background curvature for the evolution of cosmological elec- 
tromagnetic fields. The gauge-invariance of our linear equations ensures that our results are free 
from any gauge-related problems and ambiguities. We show that when the zero-order spacetime 
has open spatial sections, the magnetic flux is not always conserved. More specifically, magnetic 
fields coherent on the largest subcurvature scales are found to decay as a~ l , instead of following 
the familiar a~ 2 -law, where a is the cosmological scale factor. The reason for this deviation 
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is the general relativistic coupling between the magnetic field and the intrinsic curvature of a 
perturbed spatially open FRW universe. This magneto-geometrical interaction can change the 
evolution of the field on large scales, where curvature effects become important. The result is a 
natural superadiabatic-type amplification of cosmological magnetic fields, without the need for 
new physics and without breaking away from standard electromagnetism. 
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